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Turning the BCS-BEC crossover into a phase transition by radiation. 
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We show that the Bardeen-Cooper-Schrieffer state (BCS) and the Bose-Einstein condensation 
(BEG) sides of the BCS-BEC crossover can be rigorously distinguished from each other by the 
extrema of the spectrum of the fermionic excitations. Moreover, we demonstrate that this formal 
distinction is realized as a non-equilibrium phase transition under radio frequency radiation. The 
BEG phase remains translationally invariant, whereas the BCS phase transforms into the supersolid 
phase. For a two-dimensional system this effect occurs at arbitrary small amplitude of the radiation 
field. 
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Introduction- The Bose-Einstein condensation (BEC) 
[1] and Bardeen-Cooper-Schrieffer state (BCS) [2] are 
two extreme scenarios for the formation of the super¬ 
fluid state in fermionic systems where the only allowed 
gapless mode is the acoustic bosonic branch. In the 
BEC scenario, the fermions are first paired into compact 
two-particle complexes (molecules). These molecules ex¬ 
perience Bose-Einstein condensation with the acoustic 
low energy spectrum due to the weak repulsion between 
molecules. In the BCS scenario, weakly coupled Cooper 
pairs are formed from states near the Fermi level so 
that the characteristic size for the pair correlation signif¬ 
icantly exceeds the interparticle distance. Nevertheless, 
this weak coupling is sufficient to gap the fermionic ex¬ 
citation and leads to bosonic acoustic excitations as the 
oscillations of the order parameter. The physical effects 
occurring in between those two scenarios are referred to 
as BCS-BEC crossover. 

The BCS-BES crossover is captured by the simplest 
Hamiltonian density [3] 
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where a labels two (spin or pseudospin) states for the 
fermions described by Grassmann fields 
(summation over repeated indices is implied, and 
is the standard Pauli matrix), and the bosonic helds 
b{r), b* (r) describe the bound states of two fermions. The 
single specie energy part is described by (fi = 1) 


= (-zV - qaf / { 2 qm) + q^p - 62 , qSb, (lb) 


where the background vector a(r, t) and scalar potentials 
ip(r, t) are introduced to highlight the continuity equation 
for the total particle density n(r) = ip*tpa + 2 b* b. 

The parameter Sb describes the energy of the bound 
state when £;, > 0 or the position of the resonance when 
£(, < 0. In models of pre-formed pairs in superconductors 
[4], £{, is a material dependent parameter. In experiments 
with cold atoms Sb is the directly tunable position of 
the Feshbach resonance [5, 6]. Thus, cold atom system 


provide a versatile platform for a detailed study of the 
BEC-BCS crossover. 

In Eq. (lb), the constant A controls the coupling of the 
bound state (molecules) with fermionic continuum. If A 
is sufficiently small ^ |A|n^/^ (so-called narrow 

resonance regime), an analytic treatment of the problem 
is possible for any e^, and dimensionality D [7]. For large 
A the crossover can be investigated only numerically [8]. 

Definition of the “critical field” of the erossover - The 
arguments below are independent of the width of the res¬ 
onance. All numerical and analytical study of the ground 
state energy of the Hamiltonian (la) at :^,a = 0 indi¬ 
cates that the ground state energy density Basin, £b) is 
an analytic function of its arguments, (hence the term 
erossover rather than transition). The usual argument 
is that regardless of the values of the parameters, the 
last term in Eq. (la) leads to an anomalous average 
(i/'criT,^jcr 2 ^o- 2 ) (&) oc 6®^^. Given that apparently no 

other symmetry breakings occur, there is no sharply de¬ 
fined critical field £);(n) that separates the BEG and BGS 
regimes. Moreover, for all parameters the low energy 
excitations are described by superfluid hydrodynamics 
given by the Lagrangian 
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where dt = 9* -I- v • V is the convective derivative, 
(j) = dtO + (f and v = — a are the gauge invariant 

potential and velocities respectively, the fields n and 9 
are real, and are the Grassman fields describing 

the fermionic excitations (which for the problem of inter¬ 
est can be viewed as neutral BCS quasiparticles). This 
Lagrangian is an analytic function of variables, which ap¬ 
parently does not allow a definition of a critical field sep¬ 
arating the two regimes. The first term in Lagrangian 
(2) is protected by the gauge and Galilean invariances 
and the second term, Eas{n; Cb), is an analytic function 
of £{,. This implies that the spectrum of the bosonic ex¬ 
citations (phonons) also must be analytic. However, the 
spectrum of the fermionic excitations e(fc) experiences a 
reconstruction that allows for the definition of the critical 
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FIG. 1. The fermionic quasiparticle spectrum e(fc) from 
Eq. (2) for et close to its critical value (3a). Insets are the 
fermionic spectra deep in the BEG and BCS regimes. 


field el{n). 

In the deep BEC regime the fermions are entirely de¬ 
coupled from bosons so that the spectrum has a mini¬ 
mum at A: = 0. In the opposite limit, in the deep BCS 
regime, the spectrum of the quasiparticles has minima 
on the Fermi surface k = kp, see inset to Fig. 1. As the 
transition from point to sphere can not be analytic, there 
must exist a point el(n) such that 

dh{k-,eb,n)/dk^\^^^^ e.=e=(n) = 0 - ( 3 a) 


fermions described by the Grassmann fields f* , / which 
originally do not interact with any of the particles of the 
original problem (1). In the context of cold atom sys¬ 
tems this would be given by a third hyperfine state. The 
radiation induces transitions between the third species 
and one of the fermions from Eq. (1), so that in terms 
of the low energy theory ( 2 ) it creates or annihilates two 
fermionic excitations, 

Hrf = r + A/) / + + h.c.] , (4) 

where F{t) is proportional to the strength of the radia¬ 
tion field, and Af > 0 is the boundary of the spectrum 
of the third specie. The functional form of the second 
term in Eq. (4) is protected by gauge invariance. Let us 
concentrate on the case of the monochromatic radiation 
F{t) oc and introduce detuning as 

d = A{n,ef,) + Af — oj. (5) 

If d < 0 the Hamiltonian (4) creates real fermion pairs 
and phonons and therefore leads to entropy growth (heat¬ 
ing). For d > 0 real processes are not allowed (in fact, 
multi-photon real processes are also forbidden as the 
Hamiltonian (4) necessarily creates one / particle per one 
photon ), and the coupling (4) introduces a correction to 
the ground state energy density 


We will call el{n) from Eq. (3a) the critical field of the 
BEC-BCS transition. The fermionic spectrum for small 
momenta can be written as 

e{k, n) = A(n, ef) -I- a (e& - el) k"^ /2 F (3b) 

where a,(3 > 0. At fields below the transition Sb < el, 
the spectrum is Mexican hat shaped with the position of 
the minimum |fc| = kp, and its energy A — A: 

kp = \J{a//3) {el - e&), A = ^fc|./4 (3c) 

At first glance, the definition (3) appears to be of no 
physical consequence. Indeed, at £{, = el the fermionic 
spectrum remains gapped so that there is neither a re¬ 
construction of the ground state nor a thermodynamic 
singularity at finite temperature. However we now show 
that by an arbitrarily weak time-dependent perturbation 
it is possible to induce a spontaneous symmetry break¬ 
ing of the ground state of the two-dimensional system 
precisely at the critical point (3a). For a finite perturba¬ 
tion, the theory outlined below predicts the formation of 
the incommensurate supersolid state via a weak quantum 
first order phase transition. The periodicity of this phase 
will be determined by the “order parameter” (3c). 

Turning the crossover to the transition by coupling to 
radiation - The controlled radiative coupling to the “ex¬ 
ternal” fermions has been experimentally demonstrated, 
for example in Ref. [9]. It involves a third species of 


= m = ( 6 ) 

where ^(fc) = k^j{2m) + e{k) — A -|- d is the energy of the 
virtual state consisting of two excited fermions, and the 
meaning of the superscript ( 0 ) will become clear shortly. 

The correction SEq^ is logarithmically divergent as 
d —> 0. This corresponds to a photon with energy just 
sufficient for the excitation of / and 7 fermions with zero 
momentum. Above the BCS-BEC transition field, £& > 
el- This is the lowest energy of any excitation of a / and 
a 7 fermion and ( 6 ) is the final answer. As this energy 
correction by itself is not observable, radiation above the 
threshold, d > 0 , does not lead to any changes in the 
properties of the ground states of the system. 

The situation changes qualitatively below the BCS- 
BEC transition, £& < el- Indeed the minimal energy 
of the pair excitations is given by a fermion / with k = 0 
and one of the 7 fermions at |fc| = kp from Eq. (3c), 
i.e. the lowest boundary of the two particle continuum is 
given by A(n) -|- A^ — A. If this were to appear in the de¬ 
nominator in Eq. ( 6 ) this would give a divergence already 
at d = A, and a singular correction to the ground state. 
For the homogeneous state this is impossible as the trans¬ 
lational invariance of the ground state and of the Hamil¬ 
tonian (4) prohibits the excitation of two-quasiparticle 
with the total momentum |fc| = kp from a zero momen¬ 
tum photon. The main idea of this Letter is to propose 
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a spontaneous breaking of translational symmetry to en¬ 
able this process. 

The resulting state has no currents and the variation 
of density Sn{r) (n) is periodic in space (5n(r -|-jiti -|- 
^ 2 ^ 2 ) = <5n(r) (here ti 2 are the primitive translation vec¬ 
tors) It can therefore be classified as a supersolid state. 
If the primitive vectors of the reciprocal lattice bi ^2 have 
the length of the excitation of the lowest state be¬ 
comes allowed and the logarithmically divergent negative 
correction to the ground state is present. We will see 
that this correction can overcome the positive contribu¬ 
tion to the ground state energy from the compressibility 
{\/2){5nYE qs/ dv?)F=o thus making the supersolid 
state energetically favorable. 

Supersolid state and the phase diagram- In the pres¬ 
ence of the periodic density variation, the fermionic 
gap in Eq. (3b) also acquires a spatial variation, A = 
{dA/dn)n={n)Sn. The correction ( 6 ) changes due to the 
effect of the periodic potential produced by the super¬ 
solid: 


= - 


b,a 


keBZ 


drk 


£b,Q(k); £ 


^b,Q; (^) 


(7a) 

where the quasimomentum integration is performed 
within the first Brillouin zone, b is a vector of the recip¬ 
rocal lattice, a labels the band for the 7 fermion in the 
periodic potential, described by Schrodinger equation for 
the Bloch functions, Mj,k(r -I- t) = Mj^k(r)e**^*, 


/3(fc| + V")V4 +A(r) 






The energy of the two particle virtual state is ^b,a(k) = 
d — A -I- (k -I- b)^/ {2m) + ^ct(k), and the matrix elements 
connecting excited states to the ground state are 

F^“(k) = [ dVe-*('^+‘’)-u„,k(r). (7c) 

J UC 


The integration is within the lattice unit cell of area Sue, 
and Bloch functions are normalized as 



(7d) 



FIG. 2. The weak coupling spectrum of the 7-fermions in 
the periodic potential. The zoom shows the band structure 
near T point. The labeling corresponds to the irreducible 
representations of the point symmetry group 'Dq. Right inset: 
The basis vectors for the reciprocal lattice bi,,.,,6 and the first 
Brillouin zone. 


On symmetry grounds only Ai state, invariant under 
the symmetry group (see Fig. 2), can contribute to the 
matrix elements (7c) and 

(k = 0 ) = -V6F,A/A, U, (0) = -2A. (9) 

The linear in A shift of the lowest energy level is the 
signature of the triangular symmetry, Vq] the shift makes 
this lattice the most energetically profitable in compari¬ 
son with, e.g., square one. 

The main contribution to the energy differences be¬ 
cause the symmetry broken and the symmetric states 
comes from the lowest energy part of the spectrum. For 
the calculation with logarithmic accuracy, the partial 
contribution £ can be approximately written as £ « 
d — A — 2A -I- {2m). It yields 


^( 0 ) 


Egs = 


dA^m 

TT 



Y{d-2A) + g 


, ( 10 ) 


where the detuning from the lowest excitation energy is 
given by d = d — A, and function Y is defined as 

Y{X) = In [kl/{mX)] . (II) 


For A = 0 , equations (7) are nothing but the expres¬ 
sion ( 6 ) folded into the first Brillouin zone, as all the 
other couplings (7c) vanish. 

For small A the relevant part of the spectrum can be 
described in the almost free particle approximation. Con¬ 
sider triangular lattice 

6 

A(r) = -A ^ exp(ib, • r), (8) 

/=! 

where the vectors bi are shown on Fig. 2 a). 


The last term in Eq. (10) is the compressibility contribu¬ 
tion and the positive constant, 

g = ( 7 r/m) {d^E gs/ dn^){dA/dn)~'^ , ( 12 ) 

is of the order of unity. 

The correction given in Eq. (10) is the main result 
for the ground state energy at the lowest order in |Fp. 
It shows that the broken symmetry supersolid state is 
always energetically profitable for any finite F. How¬ 
ever the potential is apparently pathological, as Y{X) 
diverges as X —>■ 0. This infinite growth is an artifact 
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FIG. 3. The correction to the ground state energy as a func¬ 
tion of the order parameter, for liquid a) and supersolid c) 
states. Though such energy profile is typical for the first or¬ 
der phase transitions, its functional form is different from the 
Landau type expansion. 


of the lowest in F approximation, as the presence of the 
external field leads to the level repulsion. This level re¬ 
pulsion cuts the logarithm, and we turn to study of this 
effect. 

For \Fcr\ ^ A it is sufficient to take into account only 
the two lowest energy states which couple to the radia¬ 
tion and their interaction with the reference state without 
fermions. Then, the partial energy £ from Eq. (7a), with 
account of Eq. (9), becomes the lowest eigenvalues of the 
three state effective Hamiltonian 


Heff - 



0 

rf-A-2A+|^ 

%/6f;a 

A 


F A 

I , (13) 


Straightforward calculation leads to the replacement 


Y{d-2A) Y 



- 2A 


6A2|E<, 


( 10 ') 


in Eq. (10), where d{F^) = d—A-|-jF^P/A has the mean¬ 
ing of the lowest energy of the two-fermionic excitations 
shifted by the RE field. 

The resulting form of the energy profile (10), (10') is 
shown on Fig. 3. It shows two locally stable state char¬ 
acteristic of the first order phase transition. Direct in¬ 
spection shows that the supersolid state becomes more 
energetically profitable when \d{Fcr)\ < dc{Fa) where the 
critical detuning is given by 

( V3l{F,)\F,\A A^g 
[ 2A2 J ^ ^ 

The resulting phase diagram is shown on Fig. 4. 

In conclusion, we have noted the BCS-BES crossover is 
necessarily followed by the reconstruction of topology of 
the spectrum of the fermionic excitations and the criti¬ 
cal held can be rigorously dehned as the point of such 
change. We suggested an experimental scheme which 



FIG. 4. The proposed phase diagram at the fixed amplitude 
of the radiation field (main panel). The insets show the phase 
diagram at fixed position of the resonance Cb in the BEG a) 
and BGS b) regions. The “hysteresis” region of the phase 
diagram depends on the preparation of the system and may 
correspond either to super-solid or the non-stationary normal 
phase heated by r.f.-radiation. 


transmutes this reconstruction of the excitation spectrum 
into a change in the symmetry of the ground state. For 
this scenario the supersolid state is predicted to form. 

The actual process by which the supersolid state forms 
is apparently quite complex. The process is an inherently 
non-equilibrium, zero temperature and first order phase 
transition. Each of these features alone bring interest¬ 
ing facets to the issue of the phase transition kinetics. 
Therefore this transition could be an interesting arena 
for testing theories of phase transition kinetics. 
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The purpose of this supplementary section is to obtain 
explicitly express the parameters of the hydrodynamic 
description [Eqs. (2) - (3) of the main text]. We restrict 
ourselves to the two-dimensional case, = 2, in the 
narrow resonance regime, 

ri = Xy/n( — \ <C 1. (S.l) 

Vm/ 


SUPPLEMENTAL MATERIAL: 
HYDRODYNAMICS IN THE NARROW 
RESONANCE LIMIT. 


Moreover, we will consider the position of the resonance 
Eb near the critical one as discussed in the main text. 

We begin with the Hamiltonian density (1) of the main 
text, setting all the gauge fields to zero. 







(S.2) 


As we are dealing with two-dimensional systems, all 
the observable quantities expressed via the bare parame¬ 
ters of the Hamiltonian contain logarithmic divergences, 
however, the relations between different observables are 
free of such divergences. To illustrate this point, we note 
that the quantity £b is not the observable location of the 
resonance. We may calculate the physical resonance 
by computing the correction to the energy of one b par¬ 
ticle at zero momentum due to the excitations of two 
virtual fermions: 


SU (w) = -A^ 


(Pk I 

^ ' m 


Ltt \w/ 


(S.3) 


where A is an unphysical high energy cutoff. 

The physical location of the resonance is determined 
by the self-consistency equation 


mX^ 


Z I ■ (S.4) 


where A is some high-energy cut-off. The value of E,. is 
an observable position of the bound state at < 0 and 
the position of the resonance at Er > 0. 

We proceed to calculate the properties of the ground 
in the in the saddle point approximation, which is valid 
in the narrow resonance regime. We take the spatially 
homogeneous ansatz 


6 = A/AgM. (S.5) 

and introduce the thermodynamic potential density so 
that 


Basin) = H(/r) - ^a^H(/r), 


(S.6) 
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where /r(n) is found from 
n = 


(S.7) 


The the thermodynamic potential is found as the 
ground state of the mean-field version of the Hamiltonian 
(S.2) 


H = — (2/r -I- Eh) + c\ 


2m 


M Ca 


T-y J 
2 O’! a'llT2 0'2 


h.c. 


(S.8) 


After Bogoliubov rotation of the fermion operators in 
Eq. (S.8), we have their spectrum for small k and small 
fi is 


e{k) = 


2m 


|A|- 


~ k-j + A^ 


1/2 


(S.9) 


2m|A| 8m2|A| 


Comparing (S.9) with Eq. (3b) of the main text we see 
that the critical point eg is determined by /x = 0. Since 
we are investigating the vicinity of the region around this 
critical point we can restrict ourselves small /x limit. 

The thermodynamic potential at zero temperature is 
given by, 


A2 


= (-£& - 2/1) ^ 


-f 


drk 


2m 

i2 


2m 


A' 


, A^ 'ni .n. f -v/eAN 

« {-Eb - 2/i) ^ - —A In — j 

o 1 ^ ^2, f Ve\Er 

= (E. - 2/x) ^ --A Inf^ 


OifiA) 

+ 0(/iA), 
(S.IO) 


where in the last line the physical resonance Er from 
Eq. (S.4) is used to obtain the expression free of the log¬ 
arithmic divergences. 

Minimizing Eq. (S.IO) with respect to A gives, 

mX^ f \Er\\ 

j , 

with the resulting expression for the fermionic gap 
A(/x) = A(0)exp ; A(0) = |£;^|exp 

™(s.ii) 

and the thermodynamic potential 

H(/x) = -^A2(;x); (S.12) 

OTT 

Varying /x and enforcing Eq. (S.7), we obtain 

_ 2A2(^) 

^ A2 


(S.13) 


and the equation of states (S.6) as 

As noted above the critical detuning eg(n) = —Er is 
defined by /x = 0. Equations (S.ll) and (S.13) give 

which with the logarithmic accuracy yields. 


e,(n) = —l„(^ — 


(S.15) 


The values of the remaining parameters entering into 
Eq. (3b), thus are [see Eq. (S.9)] 


A ACn 1 1 

A = « = 7—r> /3 = 


^2 ’ 4mA’ Sm^A 


(S.16) 



